Update on the hadron spectrum with two flavors of staggered quarks 
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We present an update on the MILC Collaboration's light hadron spectrum calculation with two flavors of 
dynamical, staggered quarks. Results are presented for gauge couplings 5.30, 5.415, 5.50 and 5.60, with a range 
of quark masses for each value of the coupling. We present extrapolations of mjv/mp to the continuum limit for 
fixed values of m-n/mp including the physical one. 



Over the past few years the MILC Collabora- 
tion has been engaged in a series of spectrum 
calculations with dynamical quarks and in the 
quenched approximation, using both standard 
and improved actions[|]. One objective of this 
work is to obtain suflicient control over system- 
atic errors so that one can make reliable extrap- 
olations to the limits of zero lattice spacing and 
physical quark masses. Another objective is to 
make detailed comparisons between the spectra 
with quenched and dynamical quarks, and with 
standard and improved actions. Clearly, to reach 
these objectives requires high statistics calcula- 
tions for a sufficiently wide range of lattice spac- 
ings and quark masses to enable us to make ex- 
trapolations to the continuum and chiral limits. 
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It also requires that these calculations be carried 
out on large enough lattices to avoid finite size 
effects. 

In this note we provide an update on our results 
for two flavors of dynamical, staggered quarks us- 
ing the standard gauge and quark actions. Re- 
sults for the quenched approximation and with 
improved actions are presented elsewhere in these 
proceedings!^. We have carried out simulations 
at four values of the gauge coupling, 6/17^ = 5.30, 
5.415, 5.50 and 5.60 using at least four values of 
the quark mass at each coupling. Lattices were 
generated with the refreshed hybrid molecular dy- 
namics algorithm. Once the lattices were equili- 
brated the light hadron spectrum was measured 
every fifth molecular dynamics time unit. Lat- 
tices were saved after every tenth time unit for use 
in other projects. With the exeption of a few runs 
at large quark mass, measurements were made on 
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at least four hundred lattices at each values of the 
gauge couping and quark mass. 

As one would expect, the hadron masses evalu- 
ated in lattice units are strongly dependent on the 
gauge coupling and bare lattice quark mass, aniq. 
Here a is the lattice spacing. The one exception 
is the mass of the Goldstone pion which is weakly 
dependent on the gauge coupling. In Fig. 1, we 
plot the p mass as a function of the bare quark 
mass for the two weakest values of the gauge cou- 
pling that we have studied. For comparison we 
also include in this figure results from quenched 
calculations with staggered quarks at gauge cou- 
plings 5.70, 5.85 and 6.15. Notice that for fixed 
values of the gauge coupling, arup decreases more 
rapidly for decreasing aniq in full QCD than in 
the quenched approximation. This is due to the 
dependence of the lattice spacing on the quark 
mass in full QCD. 
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Figure 2. (5^ as a function of aniq for full QCD. 

have studied. The trends are as expected. The 
value of for 6/3^ = 5.60 is close to, but slightly 
below that found in the quenched approximation 
for gauge coupling 5.85. 
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Figure 1. The p mass as a function of the quark 
mass for full and quenched QCD. 

An important question regarding any simula- 
tion with staggered quarks is the extent to which 
flavor symmetry is restored as the lattice spac- 
ing and quark mass are decreased. One mea- 
sure of flavor symmetry violation is the quantity 
5^ = (m^2 — m^)/m^. In Fig. 2, we plot as 
a function of aniq for each of the couplings we 
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Figure 3. The Edinburgh plot for full and 
quenched QCD. 

In Fig. 3, we show the Edinburgh plot for our 
two weakest gauge couplings using the Goldstone 
pion in the ratio mT^/nip. Once again we include 
quenched results for gauge couplings 5.70, 5.85 
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Figure 4. mN/mp as a function of the square of 
the lattice spacing for m-^/mp = 0.5 

and 6.15 in this graph. Here too, the fuh QCD 
results at 5.60 are quite close to the quenched 
ones at 5.85. 

Because we have carried out calculations with 
a number of quark masses at each value of the 
coupling we have studied, we can perform fits to 
the hadron masses as a function of aniq in order 
to extrapolate to the chiral limit and to interpo- 
late for comparisons with quenched and improved 
action calculations. We are analyzing a variety of 
fitting functions. Here we present results from fits 
of the nucleon and rho masses to the form 

m = mo + aruq + (3mg^^ + ■yrrig. 

In Fig. 4 we show mif/mp for the four values 
of the gauge coupling we have studied, in each 
case interpolating to the value of arUq for which 
mT^/m,p = 0.5. The a:;-axis in this figure is (anip)'^ 
interpolated to the same value of the quark mass. 
This quantity gives a measure of the square of the 
lattice spacing. 

Fig. 5 is a repeat of Fig. 4, this time with ex- 
trapolations to the value of the quark mass for 
which m-^lmp takes on its experimental value, 
0.1753. The errors for the individual points in 
Figs. 4 and 5 were determined by a jackknife anal- 
ysis. The solid lines in these figures are fits to 



Figure 5. mN/mp as a function of the square 
of the lattice spacing with m-^jmp fixed at its 
experimental value. 

the four interpolated or extrapolated points of the 
form Q + /3(amp)^, where a and /? are fit parame- 
ters. The confidence levels of these fits are shown 
in the figures. The fits allow an extrapolation of 
mjq jmp to the continuum limit for fixed values of 
m-^jmp, and they provide support of the expec- 
tation that for staggered quarks the leading cor- 
rections are of order . These results are quite 
encouraging; however, a better understanding of 
the chiral extrapolation, as well as additional cal- 
culations at weaker coupling and smaller quark 
masses are needed in order to obtain definitive 
results. 

This work was supported by the Department 
of Energy and the National Science Foundation. 
Computations were carried out at ORNL, PSC, 
NCSA, lU, SDSC, CTC and MHPCC. 

REFERENCES 

1. The MILC Collaboration, C. Bernard et ai, 
Nucl. Phys. B (Proc. Suppl.) 47 (1996) 345; 
Nucl. Phys. B (Proc. Suppl.) 53 (1997) 212; 
Phys. Rev. D55 (1997) 1133. 

2. The MILC Collaboration, C. Bernard et al, 
these proceedings. 



